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ABSTRACT 

In the tidal disruption of a star by a black hole, roughly half of the stellar mass becomes bound and 
falls into the black hole, while the other half is ejected at high velocity. Several previous studies have 
considered the emission resulting from the accretion of bound material; we consider the possibility that 
the unbound debris may also radiate once it has expanded and become transparent. We show that the 
gradual energy input from hydrogen recombination compensates for adiabatic loses over significant 
expansion factors. The opacity also drops dramatically with recombination, and the internal energy 
can be radiated by means of a cooling-transparency wave propagating from the surface layers inward. 
The result is a brief optical transient occurring ~ 1 week after disruption and lasting 3-5 days with peak 
luminosities of 10 40 — 10 42 ergs s _1 , depending on the mass of the disrupted star. These recombination 
powered transients should accompany the x-ray /ultraviolet flare from the accretion of bound material, 
and so may be a useful signature for discriminating tidal disruption events, especially for lower and 
intermediate mass black holes. 
Subject headings: 



1. INTRODUCTION 

Stars unfortunate enough to pass too close to a mas- 
sive black hole are torn apart by tidal gravity forces. In 
the process, roughly half of the stellar material becomes 
bound, circularizes, and may eventually accrete onto the 
black-hole. The other half of the star is ejected from the 
system at high differential velocities. 

Tidal disruption events offers one means of probing a 
central black hole in quiescent galaxies. The dominant 
observational signature should be a bright x-ray flare 
powe red by accretion of bound material onto the black 
hole jReesfll988t lUlmerlll999t iRamirez-Ruiz fc Rosswoel 
12009( 1. Candidate disruption fla res have been de- 
tected in x-ray/ultraviolet surveys jDonlev et all 12001 
iGezari etalll2006t lEsquej et al.ll2007[[ in some case with 
counterpart emission in the optical (jGezari et al.ll2009f ). 
The light curves decline as power laws with luminosities 
comparable to the Eddington luminosity for 10 6 — 10 7 M 
black holes. Nevertheless, it can be difficult to uniquely 
identify these flares as tidal disruption events as opposed 
to other forms of nuclear activity. 

It is therefore useful to determine additional discrim- 
inating signatures of tidal disruption events. A few 
possibilities have been suggested. For deeply pene- 
trating orbits, strong tidal compression perpendicular 
to the orbital plane leads to the formation of shocks, 
which breakout out in short ( ~ 100 sec) but bright 
(L ~ 10 43 ergs s" 1 ) x-ray bursts (iKobavashi et al]l2004 : 
iBrassart &; Lum inct 2008; Guillocho n et al.ll2008l ). Later 
on, streams of material falling onto the black hole may 
undergo collisions prod ucing short UV flares of luminos- 
ity 10 40 - 10 41 ergs s" 1 (|Cannizzo et al.lll990tlKim et all 
1999). In addition, the x-ray/UV luminosity from ac- 
cretion may be absorbed and reprocessed by the un- 
bound debris into optical emis sion, leading to optical 
luminosities of ~ 10 ergs s _1 (|Bogdanovic et aLl I2004J : 
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iStrubbe fc Quataertl 120091 ). If a substantial fraction of 
the accreting material is blown in a wind, the optical ra- 
diation escaping fro m this outflow may be much brighter, 
10 42 -10 43 ergs s" 1 (jStrubbe k. Quataeru[2009( ). Optical 
emission is of particular interest, as observational surveys 
are typically most sensitive at these wavelength. 

Little attention has been given to emission from the 
unbound debris of disruption. This material, which ex- 
pands differentially with velocities ~ 10 4 km s _1 in some 
ways resembles a supernova remnant, though usually 
lacking any radioactive isotopes to power the light curve. 
While the initial internal energy of the star can be fairly 
significant (i?i n t ~ 10 48 ergs), most of it will be lost 
to adiabatic expansion before the ejecta becomes suffi- 
ciently transparent to radiate. For an ideal gas with 
adiabatic index 7 = 5/3, the energy and temperature 
evolve adiabatically as E oc p^~ x ^ oc p 2 / 3 . A decline in 
density by at least a factor of ~ 10 9 is required before 
the debris might become translucent, at which point the 
internal energy and temperature will have decreased by 
a factor ~ 10 6 . Hence the apparent lack of interest in 
diffusive thermal emission from the unbound remnants 
of tidal disruption. 

There is, however, another important energy reser- 
voir of the debris - ionization energy, which amounts to 
~ 10 46 ergs per solar mass of ionized hydrogen and is not 
subject to adiabatic degradation. When recombination 
eventually sets in at temperatures Tj ~ 10 4 K, the en- 
ergy released per atom is an order of magnitude greater 
than the mean thermal energy, Ci/kTi « 15. Recombi- 
nation is therefore a gradual process, steadily inputting 
energy to maintain the gas temperature near T^. The 
thermodynamic evolution is described by E oc p(T 3-1 ) 
where 73 is the third generalized adiabatic index, now 
a function of density and temperature. In the region of 
partial ionization 73 — 1 w 2fcT/e; w 0.1 (see 33]) so that 
the energy decline is indeed gradual. 

When the effects of ionization are included, the debris 
can expand by a much larger factor, and reach much 
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Fig. 1. — Geometrical distribution of the debris (2 days after 
pericenter passage) for a /3 = 5 tidal disruption of a solar type star 
by a 10 6 Mq black hole. The location of individual points was 
determined by following the collisionless orbits of particles under 
the gravity of the black hole, starting from a pancake configuration 
at pericenter. The top panel shows a slice through the orbital plane, 
while the bottom panel shows the distribution perpendicular to the 
orbital plane (note the vertical axis of the bottom panel is zoomed 
by a factor of 10). Black points denote bound material and red 
points unbound material. 

lower densities, before the internal energy is depleted. 
When near neutrality is reached, the opacity also drops 
sharply due to the elimination of electron scattering. It 
may be possible for radiation to escape at this time, tak- 
ing the form of a cooling-transparency wave that prop- 
agates from the surface inwards. The escape of pho- 
tons will be aided by the aspherical geometry of the de- 
bris, which has been stretched by tidal forces into a thin 
stream, with lower optical depth along the direction per- 
pendicular to the orbital plane 

It is therefore plausible that the unbound debris of tidal 
disruptions gives rise to a brief optical transient powered 
primarily by the energy tapped from recombination. We 
call such an event a recombination transient (RT). The 
light curves, even if they are dim, might be relevant for 
upcoming synoptic optical surveys which will probe for 
the first time extra-galactic transients with luminosities 
between classical novae (Mr > —10 mag) and super- 
novae (Mr < — 14 mag). Here we elaborate on the basic 
physical ideas and make approximate predictions of the 
light curves. 

2. THE DEBRIS OF TIDAL DISRUPTION 

The dynamics of tidal disruption have been 
widel y studied, both analytically and n umerically 
(e.g.. iReesI 119881: iCarter fc Luminetl 119831: iKoc hanek 
1994 iBrassart fe LuminetJ 120081: lLodato et all 120091 : 
Guillochon et alj 120081: iRamirez-Ruiz fe Rosswoel 120091 ). 
Disruption occurs when a star passes within the tidal 
radius R T = R*(M hh / M*) 1/3 , where M bh is the black 
hole mass, and M± and i?* the stellar mass and radius. 



The penetration factor (3 = Rt/R p describes the depth 
of the encounter, where R p is the radius at pericenter. 
The velocity at pericenter is v p — (2GAfbh/-R P ) 1 ^ 2 or 



6 x 10 4 (3 



x ' 2 m^W^-i/a km s" 1 (1) 



where m 6 = M bh /(10 6 M Q ), = M*/M , and r* = 
R*/ Rq . Because the initial velocity of the star at infinity 
(~ 100 km s~ ) is much less than v p , it is assumed to 
approach on an essentially parabolic orbit. 

Within Rt, the self-gravity of the star becomes sub- 
dominant. Material on the far side of the star, being 
more distant to to the black hole by possesses an 
excess specific energy relative to the center of mass, Ae = 
GMbhR+f Rp. This material becomes unbound with an 

expansion velocity at infinity Vaa — (2Ae) 1 ^ 2 or 



Voo « 6 x 10 3 (3 ml /e i 



i^^r*- 1 ' 2 km s- 1 (2) 



The half of the star nearer to the black hole, with specific 
energy — Ae, becomes bound and eventually may accrete 
onto the black hole. 

The unbound material is ejected from the system 
on hyperboli c orbits, and the debris elongates i nto a 
thin stream fe hokhlov fe Melial 119961 : iKochanekl 119941 : 
iRosswog et al.l 120091 ) . We describe the final geometry 
with the parameters R ,Rt, and R z which denote, re- 
spectively, the spatial extant of the debris in the or- 
bital direction, the transverse direction within the or- 
bital plane, and the direction perpendicular to the or- 
bital plane. The thickness of the stream within the or- 
bital plane is related to the difference in specific energy 
between the leading and trailing edges of the star as it 
approaches pericenter, Ae t = GMbhR* 2 / R p - Comparing 

the corresponding velocity at infinity v t = (2 Aet) 1 / 2 to 
Woo provides a estimate of the axis ratio E t = R /Rt- 



E t w voo/vt « 10 /r^me/m*) 1 / 6 



(3) 



In the dimension perpendicular to the orbital plane 
(the z-dire ction) tidal forces com press the star into 
a pancake (jGuillochon et al.l 120081 ). The compression 
rais es the central density b y a factor proportional to 
P 3 ([Carter fe Lumine"tl H983) , with a corresponding adia- 
batic increase in temperature of T oc p^ 7-1 ) oc (3 2 . Even- 
tually the star rebounds under the increased central pres- 
sure. For deep penetrations (f3 > 3), shocks form and 
propagate outward in the z-direction. The temperature 
of the shocked material will be of order the virial tem- 
perature 



2 GM+firrip 
R*kb 



1 x 10 7 (3 2 m+r+- 



K 



(4) 



where /i = 0.5 is the mean particle mass of ionized hy- 
drogen. The velocity in the z-direction of the shocked 
material will be of order 



v z w {ZkT+/nm p ) 1/2 w 600 (3 m 4 1/2 rr 1/2 km 



(5) 



Comparable shock vel ocities can be seen in t he hy drody- 
namical calculations of lBrassart fc LuminetJ ([2008) . This 
z- velocity serves to increase final the spatial extent of the 
unbound debris in the ^-direction. The orbital plane of 
a gas element with non-zero z- velocity at pericenter will 
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be inclined to the plane of the stellar orbit by an angle 
arcsin(ti z /w p ). The axis ratio of the debris in the direc- 
tion perpendicular to the orbital plane E z — R Q /R Z is 
then approximately 



v p /v z = 100 /T^me/m*) 1 / 3 



(6) 

We can model approximately the final shape of the 
unbound debris by simply following the orbits of colli- 
sionless test particles in the gravitational potential of the 
black hole, while ignoring self-gravity. For initial condi- 
tions at pericenter we take a uniform pancake of radius 
Rq moving in the y-direction with speed v p and with a 
z-component of velocity given by Eq. [5l Two days later, 
the final thin stream geometry (Figure [lj is roughly de- 
scribed by axis ratios of E z = 100 and E t — 10. 

The unbound debris eventually enters a homologous 
phase and the volume scales as V(t) = -R* 3 ^ 3 where we 
define the expansion parameter: 

C = (R R t R z )^±- = V -^l( Ez E t )-^ (7) 
-ft* -ft* 

We can also write ((t) = t/t e where the expansion time 
is defined 



t e = ^(E z E t )^ 

= 1120 /r 1 ™^ 176 



m* 1 /3 r+ 3/2 



E z E t 
1000 



1/3 



(8) 



Initially, the stellar mass and energy are centrally 
concentrated, but shocks largely homogenize the den- 
sity and temperature structures in the vertical direc- 
tion (iBrassart fe Luminetl 120081: iGuillochon et all 120081 : 
iRosswog et al.ll2009f ). This outward transport of entropy 
is an important factor in obtaining a bright transient. 
The total internal energy of the unbound debris after it 
has rebounded to its original volume will be of order its 
original value 



1GMV 
'2 R* 



2 x 10 4S j^mA*" 1 ergs. (9) 



where /e is a constant of order unity. By comparison, 
the total kinetic energy at infinity of the unbound debris 
is of order Eke = 1/4^**4 ~ /? 2 2 x 10 50 ergs. The 
ratio of internal to kinetic energy is initially Ei nt / Eke ~ 
10 _2 /3 -2 <C 1, which can be contrasted with that of 
shock driven Type II supernovae, -Emt/^KE ~ 1- The 
difference in initial entropy results in a distinct thermo- 
dynamic evolution and light curve for the two events. 

3. THERMODYNAMIC EVOLUTION 

Following disruption, the evolution of the expand- 
ing debris follows the thermodynamics of ionizing gase s 
(jKrishna Swamvlll96l] : iMihalas fe Weibel Mihalasfl98l . 
Initially the material is extremely optically thick and the 
evolution essentially adiabatic, dE = —PdV, where the 
pressure is 



P=(l 



NkT 1 



-aT 4 



(10) 



V 3 

with N the number of nuclei and V the volume. The 
mean ionization state is x = A+Xi, where Xi is the ion- 
ization fraction and Ai the abundance fraction of species 
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Fig. 2. — Thermodynamic evolution of a parcel of gas with initial 
temperature To = 2 X 10 7 K and density po = 1 g cm 3 . Left: 
Temperature (solid line) as a function of expansion factor in radius. 
The dashed line shows the T a (~ 2 law for a constant 7 = 5/3 
gas. The plateau in temperature starting at £ fa 10 is due to the 
energy input from recombination. Right: Evolution of the effective 
adiabatic index 73 — 1 (solid line) and the mean ionization fraction 
x (dotted line) as a function of the expansion factor. In the region 
of partial ionization, 73 — 1 « 0.1. 



i. For our purposes, only hydrogen and two ionization 
states of helium need be considered. The total internal 
energy is 



E = -(1 + x)NkT + aT 4 V + NA lXl e t 



(11) 



Excitation energy terms can be ignored. Once the tem- 
perature is high enough to excite the atoms the gas will 
quickly become ionized. 

Assuming local thermodynamic equilibrium the ioniza- 
tion state of each species is given by the Saha equation: 



1 - Xi V 



_rp3/2 



exp - 



kT 



(12) 



where Ci is a constant, and Zi is the partition function, 
which we take to be independent of temperature and 
density. 
From Eq. [Tl] we have 



(IT 



dE = NkT Y A,dx t 



-(l + x)NkT + 4aT 4 V 



kT 



(13) 



while the logarithmic derivative of Eq. rT2] gives 

dxi dx 
x 



Xi(\ Xi} 



dV 
V 



dT 
~T 



kT 



(14) 



Because the ionization energies of hydrogen and helium 
are fairly well separated, we approximate dx = dxi 

Combining Eqs.[13l[T4] and [T0l gives for adiabatic evo- 
lution 

dT , dV 

jr = (7s-l)y (15) 
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where 



73 = 1 



(l + x)(l + 4a)+E t A t f^(| + ^) 



(1 



(16) 

were a = aT 4 V/3NkT(l + x) is the ratio of radiation 
pressure to gas pressure. The sum here runs over both 
ionization states of helium, treated independently. 

When either the gas or radiation terms dominate, 
Eq. [TBI reduces as expected to 73 = 5/3 or 4/3, respec- 
tively. The terms under the sums, which incorporate the 
effects of ionization, are zero for x = and 1 but maxi- 
mal in the region of partial ionization of hydrogen. For 
Xji ~ 0.6 we get the smallest values of 73 



(73 - l)r 



khT , 
10— (1 



0.1 



T 



10 4 K 
(17) 

15 for temper a- 



where we used the fact that en/ksT 
tures T « 10 4 K. 

From Eq [17] we can estimate the amount of expansion 
the debris will undergo before recombination is complete. 
For an initial temperature of T = 10 7 K and density p = 
1 g cm 3 , gas energy dominates and one has at first T oc 
£~ 2 . The debris will expand by a factor £ » 10 before 
the temperature drops into the regime of recombination 
at T ~ 4 x 10 4 K. Hydrogen remains partially ionized for 
only a narrow temperature range, however this occurs 
over an extended period of time - to reduce T from 4 x 
10 4 K down to 1 x 10 4 requires additional expansion by 
a factor ~ 41/3(73-1).^ _ iqO. Thus, for these initial 
conditions, we anticipate total expansion of £r w 1000 
before recombination is complete. 

Figure [2] shows the thermodynamic evolution of a par- 
cel of gas with the above initial conditions, as determined 
from a direct numerical integration of Eq. 1151 The ini- 
tial T oc law plateaus at £ ~ 10, after which the 
gradual input of energy from recombination maintains 
the temperature at a near constant value. As expected, 
recombination is essentially complete by (r ~ 1000, and 
the T oc C -2 behavior resumes. Naturally, a higher ini- 
tial temperature allows for larger expansion factors, but 
so does a lower density. This is because at lower density 
recombination occurs at a lower temperature, which im- 
plies a lower 73 — 1 in the recombination regime (Eq. I17[) . 



4. OPACITY AND RADIATIVE TRANSFER 

After the debris has expanded by a factor £r ~ 1000 
and recombined, its internal energy is around 10 ergs 
per Mq of hydrogen. This is ~ 500 times greater than 
what is predicted for a constant 7 = 5/3 evolution, 
demonstrating the critical importance of the ionization 
effects. If the energy can be radiated on the time scale 
of days, a bright transient may be produced. However, 
if the debris remains optically thick the energy will lost 
quickly to expansion in the gas dominated regime. 

Radiation will preferentially escape along the z- 
direction, the shortest dimension with extant R z = 



-2/3^1/3 



The total optical depth along this di- 
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Fig. 3. — Rosseland mean opacity as a function of temperature 
for solar metallicity material with density p = 10 -10 at a time 
t = 5 x 10 sec. The thick solid line shows the total opacity, which 
is the sum of the opacity from bound-free (dashed line), electron 
scattering (thin solid line), free-free (dotted line) and line expansion 
opacity (dot-dashed). The open circles are from the tabulation of 
the Opacity Project. 




density 



t z = 7 x 10 3 krQ 2 m ir r* 



where £3 = C/10 3 an d kr is the Rosseland mean opac- 
ity in units of cm 2 g _1 . The debris is very optically 
thick for k cs w 0.4 cm 2 g _1 , the value for Thomson 
scattering in ionized hydrogen. However, the opacity 
drops dramatically with hydrogen recombination. From 
Eq. rjJl the free electron density decreases by a factor 
~ exp(— £H/2fcsT) ~ 10 -3 as the temperature changes 
from 10 4 K to 5000 K. 

Other relevant sources of opacity likewise decline 
sharply for T < 10 4 ; free-free because of the declining 
electron density and bound-free because of the reduced 
occupation number of excited levels. The H~ and molec- 
ular opacities, however, may start to become more signif- 
icant at lower temperatures. In addition, Doppler broad- 
ening in the differentially expanding debris enhances the 
line opacity from numerous iron group lines. In Figure [3] 
we plot the various sources of opacity as a function of 
temperature. The opacity reaches a minimum value of 
Kmin ~ 5 x 10" 4 cm 2 g" 1 at T w 5000 K. 

The unbound debris may thus begin to radiate effec- 
tively when the outermost layers of material have cooled 
to the "transparency" temperature T t ~ 5000 K. This 
emission should begin at roughly a time t t = (,Rt e after 
pericenter passage 



u = 13 Ca/rW 76 ™*- 1 / 3 ^ 3 / 2 



E z E t 
1000 



1/3 



days (19) 



rection is r 2 = pkrR z , or taking an average value for the 



Around this time, a sharp recombination front develops 
in the outer layers of debris - ionized material below the 
front is optically thick, while neutral material above will 
be optically thin if the densities are sufficiently low. Ra- 
diative cooling at the front promotes further recombina- 
tion, causing this surface to move progressively inward 
in mass coordinates. The so-called transparency wave 
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(|Zel'Dovich fe Raizerlfl967t IGrasberg fc Nadezhml Il976h 
allows for a more rapid release of the internal energy. 

The condition for a transparency wave to success- 
fully propagate is that the optical depth evaluated at 
the minimum opacity be of order unity. From Eq. 1181 
and K m in = 5 x 10" 4 this occurs for expansion factors 
C > 1500. The thermodynamic evolution of predicts 
expansion factors of ~ 1000 before recombination is 
complete. The coincidental near equality of these num- 
bers indicates that it is possible, though just marginally 
so, for a recombination wave to occur and radiate away 
a significant fraction of the remaining internal energy. 

The photosphere of the debris, which is coincident with 
the recombination front, cools by radiating a flux crT 4 . 
Assuming that gas interior to the front is nearly neu- 
tral at a temperature near T t , the speed at which trans- 
parency wave releases the internal energy and propagates 
inward is 



3/2nkT t 



600 Cfm* 



-i~ 3 



Tt 



5000 K 



km s 



(20) 

The timescale for the light curve to decline after peak 
will be comparable to the time it takes the transparency 
wave to propagate the distance R z 



t tw = 1.3 C 3 2 m ir r+~ 



J 100 



2/3 f E N 1/3 
10 



5000 K 



days 
(21) 



The luminosity of the light curve, as viewed perpendicu- 
lar to the orbital plane, will be roughly that of a black- 
body at T t with the given projected surface area. 



2xl0 40 C 3 V(^) (tt 



E 
100 



t 

2/3 



1/3 



Ti 



5000 K 



ergs s 
(22) 



The disruption of solar mass stars should give rise to 
rather brief, dim transients, while those of more massive 
stars will be significantly brighter given their greater i?* 
and Cr- 

5. LIGHT CURVES 

A realistic calculation of the optical emission from 
tidal disruption would require a 3-dimensional coupled 
radiation-hydrodynamics simulation. Here we perform 
simple 1-D calculations to illustrate the approximate lu- 
minosity and duration of the light curves. 

The hydrodynamic simulations of[Brassart fc Luminet] 
(2008) show that shocks tend to homogenize the debris 
structure in the z-direction. We consider models of ini- 
tially uniform temperature and density = M e j/V, 
where the ejected mass is M c - } — M*/2 and the initial 
volume V = AitR* 3 /S. The initial energy is given by 
Eq. [5] with f E — 1 or 2. To estimate the light curves for 
stars of different masses, we apply a simple mass-radius 
relation R*/R Q = (M/M Q ) - 8 . 

Radiation transport is most relevant along the z- 
direction, which we resolve with 1000 zones. The dif- 
fusion time is long in the ionized central regions, so ra- 
diation transport is only of critical importance in the 
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Fig. 4. — Bolomctric light curves of recombination transients 
from the tidal disruption of stars with masses M = 1,2,3,5,10 
and 20 Af0(from left to right, marked on figure). The calculations 
assume M b = 10 6 M Q , f) = 3, Et = 10, E z = 100 and a viewing 
angle perpendicular to the orbital plane. The left and right panels 
show calculations for two different values of the initial internal en- 
ergy (Eq. Qj). The temperature and density structures of the debris 
were assumed homogeneous. 

surface layers. We therefore calculate the adiabatic evo- 
lution as described in Sj3] while including radiative losses 
relevant for the cooling recombination wave, namely 



E, 



ad 



Tz StKpcaT 4 



(23) 



where 5t is the duration of the timestep. The isotropic 
equivalent luminosity at a given time, as seen from the 
brightest viewing angle, is 

L{t) = 4irR+ 2 C 2 E 2 J 3 E; 1/3 a J dre^T^r) (24) 

where r is the optical depth along the z-direction 
(Eq. [18|) . Given the small z-dimension, we neglect light 
crossing times. 

Figure 2] shows light curves of disruption events for 
stars of masses M = 1 — 10 M@. We have taken a pene- 
tration factor of P = 3 and axis ratios E z = 100, E t = 10. 
As expected, the light curves rise at about 3 to 5 days 
after disruption, and last around a few days. Following 
the peak there is a more slowly declining light curve tail 
in which the neutral, cold (T ~ 2000 K) debris contin- 
ues to radiate in the infrared. A proper description of 
the transport in this late phase would require inclusion 
of molecular opacities. 

A solar mass star reaches a peak luminosity of 4 x 
10 39 ergs/sec for /e = 1 and 2 x 10 40 ergs/sec for /e = 2. 
The brightnesses increases strongly with M*, mainly be- 
cause the more massive stars have larger radii and hence 
greater surface area. In addition, more massive stars 
have initially larger internal energies and lower average 
densities, which allows for larger expansion factors ($3}. 
The disruption of Af* > 10 M stars, though presum- 
ably rare events, would produce quite luminous recombi- 
nation transients, approaching the brightness of Type II 
supernovae (~ 10 42 ergs/sec) although with much shorter 
durations. 

The emission of the recombination transient will be 
anisotropic due to the asymmetry of the debris. To first 
order this angular dependence is given by the projected 
surface area of the photosphere along different viewing 
angles. For E z ~S> 1 the angular dependence is essentially 
L cx cos(i), where i is the viewing angle relative to the 
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Fig. 5. — Domain in M^-M^ space where the emission from 
the unbound debris exceeds the optical luminosity of accretion 
and fallback, as defined in the text. The recombination transient 
dominates in the region above the thick solid line (i.e., smaller 
black hole mass Mbh and larger stellar mass M*). The thin solid 
line shows the region where the time of peak for the recombination 
transient is comparable to the fallback timescale. 



direction perpendicular to the orbital plane. 

The spectrum of the recombination transient should 
be roughly blackbody at the transparency temperature 
Tf w 5000 K. The ongoing recombination may also lead 
to strong emission in H a . Thus the RT will likely be 
brightest in the R-band. The spectral features should is 
some ways resemble those of Type II supernovae, with 
P-Cygni profiles. The observed Doppler shifts, however, 
should be significantly lower, corresponding to velocities 
Vao/Ez ~ 600 (3 km s _1 . Given the extreme asymme- 
try, high levels of continuum linear polarization, near 
the limit of ~ 11% in plane parall el electron scattering 
atmospheres (jChandrasekhari fl960f ) should be observed 
from some viewing angles. 

6. PROSPECTS FOR DETECTION 
6.1. Domain of Astrophysical Relevance 

The recombination transient is just one of several 
sources of electromagnetic radiation in a tidal disruption 
event. An important question is whether the luminosity 
of the RT can compete with the other emission mech- 
anisms, in particular the fallback of the bound stellar 
debris and its subsequent accretion onto the black hole. 

The energy radiated during the fallback phase may 
dominate the luminosity at early times. The most tightly 
bound debris moves on an elliptical orbit with semi- 
major axis a = ^R p (R p / R+) and returns to pericenter 

on the orbital timescale £fb = 27r(a 3 /GMbh) 1 ^ 2 , or 



tfb « 40 (3 



-3 1/2 
TOfi r. 



3 / 2 m" 1 



days, 



(25) 



with less bound material falling back continually there- 
after. The streams of fallback material undergo collisions 
which, after a few orbital timescales, circularize the or- 
bits at a radius ~ 2R p dete rmined by total angular mo - 
mcntum conservation {R amircz- Ruiz fc Ross"wol[2009| ). 
The final binding energy of this circular orbit is much 
less than the binding energy of the initially elliptical or- 
bits. Assuming the energy difference is radiated with 
efficiency e over a few orbital timescales, ~ 3ifb, the lu- 



minosity is would be 

1 GM hh M+ 1 

2 2R„ 2 3tfb 



L 



fb 



4.6 x 10 43 /3 4 eo.iml /6 m+ 7/3 r*- 5/2 



(26) 



erg s 



where eo.i = O.le is highly uncertain. Assuming the pho- 
tosphere is near 2R p , the emission temperature is 



Tfb ~ 2.4 x 10 5 (3 3/2 m. 



-1/8^ 3/4 - 



9/16 



K. 



(27) 

It is possible that some of the energ y available in fallback 
goes into driving a mass outflow (|Strubbe &: Quataertl 
2009). This matter adiabatically expands to large radii 
until becoming optically thin. For certain choices of the 
fraction of mass ejected, the resulting optical luminosity 
can be quite bright L ~ 10 42 — 10 43 ergs s~ 4 . 

After circularization, the fallback material forms a 
torus which will spread out viscously and accrete onto the 
black hole. The timescale for this stage of the evolution 
is de t ermined by the viscous timescale (|Cannizzo et al.l 
fl^luln^[i999h 



7:51 or 



h 

OA 



-^rj/ 2 days. 



(28) 



where a is the stan dard viscosity parameter 
(jShakura fc Sunvaevj |l973), and h is the ratio of 
the disk height to radius and is approximately of order 
unity for thick disks. The emergent radiation may have 
a strong dependence on the structure and orientation of 
the disk. As a simple estimate we consider the luminos- 
ity to be at the Eddington limit L ct \ ft 
and emanating from the radius 2R p . 
accretion temperature 



10 to 6 ergs s " 
This gives an 



1.6 x 10 5 /^m^W/V- 1 / 2 K. 



(29) 



Both the fallback and the accretion temperature depend 
weakly on Afbh- In the case where emission peaks in the 
soft x-rays, we will observe in the optical the Raleigh- 
Jeans tail of the blackbody, where the flux is reduced by 
a factor {T^/T t )~ 2 from its peak value. 

From the above discussion, we can estimate the con- 
ditions under which the recombination transient is a rel- 
evant source of emission. Using Eq. [24] and assuming a 
simple mass-radius relation r± — m* 0,8 , we find that the 
optical luminosity of the RT exceeds that of fallback and 
accretion, respectively, when the mass of the disrupted 
star exceeds 



m* > 1.8 (/J/S) 15 / 28 ^ 5 / 28 ^ 57112 
for L t > L fb (T fb /5000 K)" 2 



and 



to* 



> 



1.4 (/3/3)- 15 / 17 to 25/34 



for L t > L a cc(T acc /5000 K)" 



(30) 



(31) 



Figure [5] depicts the relevant M*-Mh domain where the 
RT is likely to contribute significantly to the optical lu- 
minosity. While these estimates are subject to several 
uncertainties, it appears that the RT of to* s=s 1 dis- 
ruptions will be most relevant when Mbh < 10 6 M Q . 
For more massive black holes, disruption of larger stars 
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may be necessary. The RT appears to be compa- 
rable to or brighter than the optical emission result- 
ing from the repr ocessing of x-rays in the unbound 
debris, studied by iStrubbe fc Quataertl (|2009l ). How- 
ever, if the super-Eddingtq n outflows also discussed by 
IStrubbe fc Quataertl (2009) occur, their optical lumi- 
nosity can, for certain wind parameters, be as high as 
~ 10 42 — 10 43 ergs, which would completely dominate 
the RT of even massive stars, at least at early times. 

Even if the RT luminosity falls below the other sources, 
it may be discernible if it is separated in time. From 
Eq. [TH the ratio of the RT timescale, tt to the accre- 
tion/fallback timescales, ~ 3ifb, is 



3ifb 



0.1 /? 2 TO~ 2/3 m* 2/3 . 



there is solid kinematical evidence. The galaxies in ques- 
tion are dwarf elliptical (dE) galaxies, which are the most 
numerous type of galaxy in the Universe. Tidal disrup- 
tion events are also more common in dE galaxies, and 
if they contain nuclear black holes they should dominate 
the total tidal flaring rate. Assuming th at only nucleated 
dE ga laxies (dEn) contain black holes, iWang fc Merrittl 
(2004) found a total tidal disruption for dEn in Virgo 
of ~ 0.2 yr" 1 . The overall rates in the Coma Cluster 
are expected to be significantly larger. These event rates 
will be even higher if moderately massive black holes are 
present in every dE galaxy or in bulges of late-type spi- 
rals. 

Suggestive evidence has accumulated that intermediate 
(32) mass black holes (Mb h = 10 3 — 10^ Mq,) exist in some 



Thus for black holes M bh > 10 6 M Q the RT might be 
observable as a precursor to the black hole accretion and 
fallback signatures (Figure [5]). In this case it may serve 
as a confirming signature that the event being witnessed 
is in fact a tidal disruption event. 

The detection of disruption transients also requires 
that the optical emission reach some significant fraction 
of the luminosity of the host galaxy nucleus. At a dis- 
tance of 400 Mpc, ground based optical surveys with a 
resoluti on of a ~ 1" will only res olve a kiloparsec-sized 
bulge. iMcLure fc Dunlopl (|2002l) give a empirical cor- 
relation between bulge luminosity and central black 
hole mass. For Mbh = 10 s Mq, the bulge luminosity 
Lb ss 2 x 10 9 Lq, and only the disruption of > 10 M Q 
stars could contribute > 10% of the core brightness for 
galaxies. For Mbh ^ 10 7 M©, on the other hand, the 
disruption of > 3 stars would outshine the galaxy 
and be easily detectable. However, firm evidence for the 
assumed Mbh — correlation below Mbh < 10 7 M3 ls 
still lacking. 

6.2. Detection Rates 

Observational surveys are being designed which will 
probe the magnitude range (Mboi = —11 to —14) of the 
recombination transients. LSST, for example, will reach 
an R-band magnitude of Mr = 24 mag per shot. In this 
case, a M = 1 M Q disruption would be detectable out 
to ~ 100 Mpc, the distance to the Coma cluster, while 
an M = 10 Mq disruption would be detectable much 
f urther out, to z ~ 0.1. 

IWang fe Merritti (|2004l ) derived a total disruption rate 
for early-type galaxies and bulges by combining an ex- 
pression for the disruption rate 



N « 2.2 x 10" 4 yr- 1 h-°- 25 (L g /W 10 h- 2 L @ )- - 295 (33) 

with lFerguson fc Sandage! ([1991) E+SO luminosity func- 
tion for the galaxy luminosities L g . They found a rate per 
unit volume of ~ 10" 5 yr -1 Mpc" 3 , with galaxies with 
L g < 10 9 Lq (M bh < 10 8 M o ) dominating the consump- 
tion rate. Assuming that galaxies with L g -C 10 9 Lq 
hosting black holes are at least as common, we estimate 
that LSST may detect RTs of more massive stars at a 
rate of > 3OO'0>5 yr -1 , where ip > 5 is the fraction of total 
disruptions of stars with to* > 5. 

The major uncertainty in this rate is whether low lu- 
minosity galaxies contain a nuclear black hole, as M32 is 
currently the faintest system (L g ~ 4 x 1Q & Lq) for which 



globular clusters (e.g. iNovola et al.ll2008l ). The rate of 
tidal disruption events in these syste ms is expected to 
be ~ 10~ 7 yr -1 per globular cluster (|Baumgardt et al.l 
2004) . Taking a globular cluster space density of n gc ~ 

4 Mpc" 3 (<Brodie fc Straderll2006T ) we estimate the den- 
sity rate of RT to be at most ~4x lO^yr" 1 Mpc" 3 . 
In this case, LSST should detect RT events arising from 
the disruption of a 1 Mq star with a rate of ~ 5 , 0>iyr~ 1 . 
The globular clu ster luminosity function in the galaxies 
studied to date (|Brodie fc Straderl 1200(1 is well fit by 
a Gaussian distribution with a peak My ~ —7.4 and 
a ~ 1.4. The RT should thus easily outshine the clus- 
ter's optical luminosity. Non-detection of flares after a 
few years would argue against the existence of interme- 
diate mass black holes in globular clusters. 

7. SUMMARY AND CONCLUSIONS 

The purpose of this paper is has been to explore 
whether the unbound debris of tidal disruption might 
give rise to luminous optical emission, and to determine 
under what conditions this emission may be observation- 
ally relevant. The analytical results predict a brief (3-5 
day) optical transient with peak luminosities in the range 
2 x 1O 4O (M*/M ) 5 / 2 ergs s" 1 . The energy for the emis- 
sion is stored mostly in the ionization energy of hydro- 
gen, and is released when the ejecta becomes neutral and 
transparent at temperatures T « 5000 K. This physics 
will be relevant not only for the tidal disruption of stars 
by massive black holes, but to the weak ejection of hy- 
drogen envelopes in general by e.g., stellar collisions or 
pulsations. 

The predicted bolometric luminosities of the recombi- 
nation transients fall into the interesting - that is to say, 
relatively unexplored - range between novae and super- 
novae. Upcoming synoptic surveys will begin to probe 
these sorts of intermediate events. The RT of 1 Mq 
stars are relatively dim events which would only be visi- 
ble in the nearby universe (~ 100 Mpc). The disruption 
of more massive stars (M > 5 Mq), on the other hand, 
can be quite bright. Such events may not be as rare as 
one might first expect. There is strong evidence that the 
initial mass f unction of stars w ithin the Galactic center 
is top heavy ([Alexa nder 20051). Perhaps massive stars 
are more common in the cores of other galaxies as well. 

For massive black holes, Mbh > 10 6 M Q , the optical 
luminosity from the accretion of bound material likely 
surpasses the emission from the unbound debris. Never- 
theless, the RT may still be discernible as a precursor, 
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and may provide a important complementary signature. 
While the properties of the accretion flare depends most 
sensitively on the black hole mass, the luminosity of the 
RT depends primarily on the properties of the disrupted 
star, with little dependence on Mbh- In principle obser- 
vations of both components would allow one to solve for 
the parameters of the event, providing insight into both 
the black hole mass and stellar population in the cen- 
ters of distant galaxies. More basically, detection of a 
RT would provide confirmation that the event witnessed 
was in fact a legitimate tidal disruption, 

Recombination transients appear to be of greatest in- 
terest for probing lower mass black holes (Mh < 1O 6 M ). 
In these systems, the optical luminosity of the transient is 
likely to exceed that of accretion, and should outshine the 
nucleus of the galaxy. This offers one means of searching 
for low mass black holes in dwarf galaxies, or for interme- 
diate mass black hole in globular clusters. If black holes 



do lurk in these numerous locations, the recombination 
transient might be the most prominent optical signature 
of a stellar disruption. 
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